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ABSTRACT 

The WD+MS channel of the single-degenerate scenario is currently favourable 
for progenitors of type la supernovae (SNe la). Incorporating the results of 
detailed binary evolution calculations for this channel into the latest version of a 
binary population synthesis code, I obtained the distributions of many properties 
of the companion stars at the moment of SN explosion. The properties can be 
verified by future observations. 
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1. Introduction 

Type la supernovae (SNe la) have been used as a calibrated candle to probe the dynam- 
ics of the universe, leading to a significant progres s in cosmology, i.e. the determination of 
A and Q (IRiess et al.lll998t iPerlmutter et al.lll999l ). They are believed to be thermonuclear 
explosions of carbon-oxygen (CO) white dwarfs (WDs). However the nature of their progen- 
itors has remained unclear, and this still raises doubts as to the calibration which is purely 
empirical and based on the nearby SN la sample. Based on the characteristics of observed 
SNe la (e.g. light curves, chemical stratification), it seems most likely that they occur when 
the accreting CO WDs reach the Chandrasekhar limit, and sub-Chandrasekhar models ap- 
pear not to be consistent with the observations of SNe la. Chandrasehar-mass WDs can 
be created through the singl e-degenerate scenario, where the CO WD accretes mass from a 
non-d egenerate companion (INomoto. Thielemann fc Yokoil Il984i lHachisu. Kato fc Nomoto 
19991 ). or the double degenerate scenario, where two CO WDs with a total mass larger than 
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the Chandrasekhar mass coalesce Jlben fc Tutuk ovl IliiS Iwe bbink fc Ibenl [1987^. Recent 
observations indicate that there is a wide spread of delay time between star formation and 
SN la explosion, and this imply that there exist two populations of progenitors, a "prompt" 
one with a delay time less than ~ 0-1 G yr and a "tardy" one with a delay time of ~ 3 Gyr 
jMannucci. Delia Valle fc Panagialbood ). 



Usin g the archival data prio r to explosion from Hubble Space Telescope (HST) and 
Chandra, IVoss fc Nelemand (120081 ) tried to search for a possible progenitor of SN2007on (type 
la) in elliptical galaxy NGC 1404. They discovered an object at the position of SN2007on. 
The X-ray luminosity of the object and the non-detection in the optical images are fully 
consistent with the single-degenerate scenario, and the the discovery t herefore favours the 
single-degenerate scenario. However, new Chandra X-ray observations (IRoelofs et al.ll2008l ) 
and detailed astrometry of the site of SN2007on showed that there appears to be an offset 
between the SN and the X-ray source, and the offset means a small probability (~ 1%) of 
the X-ray source being related to the SN. However, the X-ray source is not unconnected 
with the SN, as the X-ray source has dimmed after the explosion and the source before 
the explosion showed an excess of soft X-ray photons relativ e to the other s o urces in the 
field. Kepler's 1604 supernova is suggested to be of type la. iReynolds et al.l (120071 ) made 
a deep Chandra observation of Kepler's supernova remnant and concluded that Kepler's 
SN is a SN la with circumstellar me dium (CSM) intera c tion. This might indicate that the 
progenitor is a massive (young) star. iMaoz fc Mannuccil (120081 ) searched HST pre-explosion 
images for NGC 1316, a radio galaxy in the Fornax cluster and a prolific producer of SNe 
la, for evolved intermediate-mass progenitor stars. The pre-explosion images (3 years before 
explosion) of the sites of SN2006dd (a normal SN la) and SN2006mr (likely a subluminous SN 
la) show no potential luminous stellar progenitors. More effort is still needed in identifying 
the progenitors of SNe la. 

Among various possible progenitor models, the WD+MS channel of the single-degenerate 
scenario is the most widely accepted one, in which a CO WD accretes mass from its 
main-sequence (MS) star or a slightly evolved subgia nt in a close binary system until it 



reaches a mass of ~ 1.378 Mm and explodes as a SN la (jNomoto. T hiele mann &: Y okoi 



1984 



Li fc van den Heuvellll997i : lHachisu. Kato fc Nomotolll99 & lLanger et aU2000l ; lHan fc Podsiadlowskil 



2004J). The companion star should survive the explosion and show distinguishing properties, 
and it is therefore a promising method to test progenitor models by identifying the surviving 
companion stars of SNe la. 



1 Note, it is quite likely that the merger product experiences core collapse rather than a thermonuclear 
explosion. 



- 3- 



Tycho Brahe's 1572 supernova is a Galactic SN la. iRuiz-Lapuente et al.l (12004 ) found in 
the remnant region that Tycho G, a star similar to the Sun but with a lower gravity, moves at 
more than three times the mean velocity of the stars there. They argued that Tycho G could 
be the surviving companion of the supernova^. Indeed, a surviving companion would have 
a high space velocity and evolve to a WD finally, and the single-degenerate scenario could 
potentially explain the proper ties of halo WDs observed by lOppenheimer et al.l (120011). e.g. 



their space density and ages ([Hansen! 120031 ; iBergeronl 120031) . Furthermore, lJustham et al. 



( 120081 ) argued that the ultra-cool WDs observed by IWolj (120051 ) might have been formed via 
the single-degenerate scenario, i.e. they might be the remnants of the non-degenerate donor 
stars. Note, however, there have been no conclusive proof yet that any individual object is 
the surviving companion of a SN la. 



Langer et al.l (120001 ) did detailed binary evolution calculations for part of the WD+MS 
channel of the single-degenerate scenario, and presented the properties, e.g. orbital velocities, 
luminosities, effective temperatures, of the companion stars at the moment of SN la explosion 
for the sample WD binaries they studied, in which the Roche lobe overflow (RLOF) starts 
when the comp anion star is in MS. However, the distri butions of the properties have not 
been obtained. ICanal. Mendez fc Ruiz-Lapuentd (120011 ) run a Monte Carlo scenario code 
to calculate the distributions of masses, luminosities and velocities of the companion stars 
of SNe la. However, they have not done detailed binary e yolution calculations for th e 
production of SNe la, and this results in big uncertainties (see lHan fc Podsiadlowskill2004l ). 
e.g. the orbital velocities obtained for the companion stars from the WD+MS channel are 
over 450 km/s, which is far too high. 



Han fc Podsiadlowskil (120041 . hereafter HP04) carried out detailed binary evolution cal- 



culations for the WD+MS channel for about 2300 close WD binaries, in which RLOF starts 
when the companion star is in MS or in Hertzsprung gap (i.e. slightly evolved). The study is 
comprehensive, and various properties of the companion stars were obtained but not sorted 
for publishing. In this Letter, I extract the properties from the data files of the calculations 
and incorporate them into the latest version of the binary population synthesis (BPS) code 



deve 



1995 



oped for the study o 



Han et al. 



1995 



Han fc Podsiadlowski 



various binary-related obje cts (IHan. Podsiadlowski Eggleton 



ianlll998l: lHan et al.ll2002l . l2003h . including the progenitors of SNe la 



2004 120061 ). and obtain the distributions of the properties. 



2 However, Fuhrmann (2005) argued that there exists a possibility of Tycho G being a thick disk star 
coincidentally passing in the vicinity of the remnant of SN 1572. 
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2. The distributions of properties of the companion stars 

In the single-degenerate scenario, the progenitor of a SN la is a close WD binary syste m 



which has most likely emerged from common envelope (CE) evolution (jPaczyriskil Il976l ) of 
a giant binary system. During the CE evolution, the envelope engulfs the core (here a CO 
WD) of the giant and the secondary, and the orbital energy released in the spiral-in process 
(i.e. orbital decay) is used to overcome the binding energy of the CE in order to eject 
it. For the CE evolution, I have two parameters: «ce the CE ejection efficiency, i.e. the 
fraction of the released orbital energy used to overcome the binding energy, and a t h, which 
defines the fraction of the thermal energy contributing to the binding energy of the CE. As 
in previous studies, I adopted «ce = «th = 1-0, which gives good matches between theory 
and observations for many binary-related objects [§. 

To obtain the distributions of properties of companion stars at the moment of SN 
explosion, I have performed a detailed Monte Carlo simulation with the latest version of the 
BPS code. The code follows the evolution of binaries with their properties being recorded 
at every step. If a binary system evolves to a WD+MS system, and if the system, at the 
beginning of the RLOF phase, is located in the SN la production regions in the plane of 
(logP 1 , Mj) for its M^, where P 1 , M\ and M^ D are, respectively, the orbital period, the 
secondary's mass and the WD's mass of the WD+MS system at the beginning of the RLOF 
(see Fig. 3 of HP04), I assume that a SN la is resulted, and the properties of the WD 
binary at the moment of SN explosion are obtained by interpolation in the 3-dimensional 
grid (M^ D , M|, logP 1 ) of the ~ 2300 close WD binaries calculated in HP04. 

In the simulation, I follow the evolution of 100 million sample binaries according to grids 
of stellar models of metallicity Z = 0.02 and the evolution channels leading to SNe la as de- 



scrib ed in HP04. I adopted the following input for the simulation (see lHan. Podsiadlowski fc Eggleton 



19951 ). (1) The star-formation rat e (SFR) is take n to b e constant over the last 15Gyr. (2) 



The initial mass function (IMF) of iMiller fc Scald (119791 ) is adopted. (3) The mass-ratio dis 



tribution is taken to be constant. (4) The distribution of separations is taken to be constant 
in log a for wide binaries, where a is the orbital separation. (5) The orbits are assumed to 
be circular. 

The simulation gives current-epoch-snapshot distributions of many properties of com- 
panion stars at the moment of SN explosion, e.g. the masses, the orbital periods, the orbital 
separations, the orbital velocities, the effective temperatures, the luminosities, the surface 



The prescription adopted here for the CE evolution is diffe r ent fr o m the A presc ri ption , but 



appears to be more physical. See 



Podsiadlowski. Rappaport k. Hani ( 2003h for details 



tor tne ufj evolution is diner ent rr o m tne A pre scription, Dut 
Han. Podsiadlowski fc Eggleton! (jl995h . bewi fc Taurisl (|200C( ) and 
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gravities, the surface abundances, the mass transfer rates, the mass loss rates of the opti- 
cally thick stellar winds. The simulation also shows the initial parameters of the primordial 
binaries and the WD binaries that lead to SNe la. Figs [T]- [7] are selected distributions that 
may be helpful to identify the progenitors of SNe la. 



Discussion 



Figs [T] and [2] are the distributions of the masses, the orbital velocities @, the effective 
temperatures and the surface gravities of companion stars at the moment of SN explosion. 
Tycho G was taken a s the surviving companion star of Tycho Brahe's 1572 supernova by 
Ruiz-Lapuente et al.l (120041 ) . It has a space velocity of 136km/s, more than 3 times the 
mean velocity there. Its surface gravity is log (g/cms -2 ) = 3.5 ± 0.5, while the effective 
temperature is T e ff = 5750 ± 250K. The parameters are compatible with Figs [1] and [2j 
As from our simulation, the recorded properties at each step show that a primordial binary 
system with a primary mass M\\ ~ 4 — 5.5M , a secondary mass M 2 \ ~ 2 — 3M Q and an 
orbital period Pi ~ 100 — 250 d would evolve to a close WD binary system with a WD mass 
Myy D ~ 0.8— 1.2M Q , a secondary mass M\ ~ 2— 3M , and an orbital period P 1 ~ 1—4 d. The 
WD binary results in SN la explosion with companion parameters (T e g and logg actually) 
in the range of Tycho G. 

However, Figs [1] and [2] are for tha t at the moment of SN explos i on, t he distributions 
could be modified due to the explosion. iMarietta. Burrows fc Fryxelll (120001 ) presented sev- 
eral high-resolution two-dimensional numerical simulations of the impacts of SN la explosion 
with companions. The impact make the companio n in the WD+MS ch a nnel lose a mass of 



0.15-0.17M©, and receive a kick of 49 - 86km/s. iMene. Chen & Hani (120071 ) adopted the 



simple analytic method of Wheeler. Lecar fc McKeel ( 19751 ). and calculated the impact to 
survey the influence of the initial parameters of the progenitor's systems. With detailed stel- 
lar models and realistic separations that were obtained from binary evolution, they obtained 
an even lower 'stripped mass', 0.03-0.13M Q , but a similar kick velocity 30-90km/s, which is 
perpendicular to the orbital velocity. A surviving companion star therefore has a mass lower 
by ~ 0.1M© and a space velocity larger by ~ 10% than that in Fig. [TJ 

The companion stars are out of thermal equilibrium at the moment of SN explosion. 
The equilibrium radii are typically larger by ~ 50% than that at the moment of SN ex- 
plosion. Therefore the surface gravity at equilibrium should be lower than that in Fig. [2j 



4 The Chandrasekhar-mass WD has an orbital velocity of ~ 50 to ~ 200km/s for a corresponding com- 
panion star's mass of ~ 0.6 to ~ 2.0M© at the moment of SN explosion. 
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Podsiadlowskil (120031 ) systematically explored the evolution and appearance of a typical com- 
panion star that has been stripped and heated by the supernova interaction during the 
post-impact re-equilibrium phase. Such a star may be significantly overluminous or under- 
luminous. Fig. [2] could be a starting point for further studies of this kind. 

Fig. [3] is the distribution of orbital periods of the WD+MS systems at the moment of 
SN explosion. If I assume that the companion stars co-rotate with their orbits, I obtain 
their distributions of equatorial rotational velocities (see Fig. 0J. We see that the surviving 
companion stars are fast rotators and their spectral lines should be broadened noticeably. 

Fig. [5]shows the distribution of masses lost during opti cally thick stell a r win d phase. We 
see that a significant amount of mass is lost in the wind. iBadenes et al.l (120071 ) found that 
the wind with a velocity above 200km/s, which is believed to be reasonable, would excavate 
a large low-density cavity around its progenitor. However, the fundamental properties of 
seven young SN la re mnants (Kepler , Tych o, SN 1006, 0509-67.5, 0519-69.0, N103B and 
SN 1885) obtained by IBadenes et al.l (120071 ) are incompatible with such large cavities. A 
l ower wind velocity or the consideration of WD rotation could help to solve the problem 
(IBadenes et alil2007f ). 



Fig. M presents the distribution of mass transfer rates at the moment of SN explosion. 
The mass transfer rates can be converted to X-ray luminosities of the systems by Lx ~ e|M|, 
where e = 7 x 10 18 erg/g gives the approximate amount of energy obtained per gram of 
hydrogen burnt into helium or carbon/oxygen. The luminosity of the X-ray source close to 
the site of SN 2007on (4 y ears before the explosion) was estimated to be (3.3±1.5) x 10 37 erg/s 
(jVoss fc Nelemansll2008l ). corresponding to a mass accretion rate of ~ lO~ 7 M /yr, which is 
consistent with Fig. [6l 

QU Carinae, a cataclysmic variable (CV), is suspected to be a SN la progenitor. It has 
a WD mass of ~ 1.2M , a once-reported orbit al period of 0.45 d, and more import antly, a 
very high mass transfer rate of ~ lO~ 7 M /yr (jKafka. Anderson fc Honeycuttl 120081 ). These 
properties are consistent with Figs [3] and [61 However, BF Eridani, ano ther CV with Mwd ~ 



1.28M , M 2 ~ O.52M and P ~ 0.27 d (INeustroev fc Zharikovl 120081 ). appears not to be a 
SN la progenitor. 

Fig. [7] is the distribution of the surface nitrogen mass fraction of companion stars at 
the moment of explosion. We see that nitrogen can be significantly overabundant (with 
corresponding underabundance of carbon due to the CN-cycle). However, the surface can 
be seriously polluted by the ejecta of SN explosions. 



The simulation in this Letter was made with acE = cuth = 1.0. If I adopt a lower value 
for a t h, say, 0.1, the birth rate of SNe la would be higher (a factor of 1.7) and the delay time 
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from the star formation to SN explosion is shorter. This is because that binaries resulted 
from CE ejections tend to have shorter orbital periods for a small a t h and are more likely to 
locate in the SN la production region (see Fig. 3 of HP04). The companion stars with orbital 
velocity V < llOkm/s would be absent from Fig. [U and the ones with log(g/cm s -2 ) < 3.1 
would be absent from Fig. [21 This is due to that WD binaries with long orbital periods are 
absent due to a small a t h- 

The distributions are snapshots at current epoch for a constant SFR. For a single star 
burst, most of the SN explosions occur between 0.1 and 1 Gyr after the burst (see Fig. 7 
of HP04). The evolution of progenitor properties with time can be understood via Fig. 3 
of HP04. A delay time from 0.1 to 1 Gyr corresponds to M\ of ~ 3.2M to ~ 1.8M , and 
to M^ D of ~ 1.2M Q to ~ O.67M for the WD+MS system, respectively. As seen from Fig. 
3 of HP04, the range of the orbital periods becomes narrower from a delay time of 0.1 to 
1.0 Gyr. Those WD binaries result in SN la explosions via RLOF. Consequently, the range of 
progenitor properties at the moment of SN la, e.g. the orbital velocities of companion stars, 
the surface gravities, the equatorial rotational velocities, the mass transfer rates, becomes 
narrower with time. The mass transfer rate would be smaller with time as M\ becomes 
smaller. 
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Fig. 1. — A snapshot distribution of companion stars in the plane of (V, M| N ) at current 
epoch, where V is the orbital velocity and M| N the mass at the moment of SN explosion. 
The number density decreases from inner regions to outer regions. Regions, from inside to 
outside with corresponding gradational grey or color in the legend (from bottom to top), 
together with the regions with higher number densities contain 50.0%, 68.3%, 95.4%, and 
99.7% of all the systems, respectively. 




Fig. 2. — Similar to Fig. [TJ but in the plane of (logT eff , logg), where T e g is the effective 
temperature of companion stars at the momen t of SN explosion, log q th e surface gravity. 
The error bars denote the location of Tycho G (iRuiz-Lapuente et al.l I2004T ) . 
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Fig. 3. — Similar to Fig. [T], but in the plane of (logP SN , M| N ), where P SN is the orbital 
period at the moment of SN explosion. 




Fig. 4. — Similar to Fig. (TJ but in the plane of (V T0 t, ^f| N ), where V vot is the equatorial 
rotational velocity of companion stars at the moment of SN explosion. 
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Fig. 5. — Similar to Fig. CD, but in the plane of (AM win d, Mf N ), where AM win d is the mass 
lost during the optically thick stellar wind phase of the WD binaries. 




Fig. 6. — Similar to Fig. [T], but in the plane of (log(— M2), M^), where — M2 is the mass 
transfer rate at the moment of SN explosion. 
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Fig. 7. — Similar to Fig. [p, but in the plane of (X N , M| N ), where X N is the nitrogen mass 
fraction at the surface of companion stars at the moment of SN explosion. 



